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INCLUSION OF EXPLICIT THERMAL REQUIREMENTS 
IN OPTIMUM STRUCTURAL DESIGN 


By Howard M. Adelman and Patricia L. Sawyer 
Langley Research Center 

SUMMARY 

A finite-element based procedure is described for obtaining minimum mass 
designs of structures subjected to combined thermal and mechanical loading and 
both strength and thermal constraints. The procedure is based on a mathemati- 
cal programming method using the Sequence of Unconstrained Minimizations Tech- 
nique (SUMT) in which design requirements are incorporated by an exterior 
penalty function. Temperatures are treated as behavior variables rather than 
fixed load-type quantities and are recalculated during the resizing process 
using a finite element thermal analysis. The procedure is limited to steady- 
state temperatures which are controlled by structural sizing only. The optimi 
zation procedure is demonstrated by the design of a structural wing box with 
both mechanical loading and external heating, subject to design constraints on 
stress, minimum gage, and temperature. The final design for these conditions 
is compared with a corresponding design in which temperature constraints are 
omitted. Temperature constraints have a significant effect on both the distri 
bution of structural material and the total mass in the final design. 

Some additional developments beyond the scope of the present work but 
needed for design of practical aerospace structures under realistic load 
situations are identified. Among these developments are design-oriented 
transient thermal analysis capability, treatment of time-dependent constraints 
and incorporation of effective temperature control devices and thermal -stress 
relieving mechanisms. 


INTRODUCTION 

The necessity for light-weight aerospace vehicles to withstand both 
mechanical loading and severe thermal environments during atmospheric entry 



and hypersonic cruise has led to the problem of designing minimum-mass struc- 
tures that must satisfy both mechanical and thermal design requirements. 
Mechanical requirements to Insure structural Integrity include limits on stress, 
buckling loads and flutter speed. Thermal requirements to avoid material 
property degradation and excessive temperatures in cargo and crew compartments 
consist of upper limits on structural temperatures. 

Previous structural synthesis work has generally concentrated on designing 
structures under mechanical loads for mechanical requirements. Thermal loads 
were incorporated as specified structural temperatures supplied by an indepen- 
dent thermal analysis and equivalent loads were added to the applied mechanical 
loads (ref. 1-3). Temperatures entered into the synthesis problem indirectly 
through the strength requirements and explicit thermal requirements were not 
included. The temperatures of points on the structures were not recalculated 
as the structure was resized. Thus, good initial estimates of the structural 
member sizes were required to avoid significant differences between the tem- 
perature distributions in the initial and final designs. 

Two logical improvements to existing thermal -structural design methods are 
incorporation of temperature constraints and recalculation of temperatures dur- 
ing resizing. Early recognition of the need for these features is reflected in 
reference 4 which describes the optimization of an ablating heat shield with 
upper limits on temperature, and reference 5 which indicates the importance of 
recalculating temperatures of structural elements during resizing. The purpose 
of the present work is to combine the aforementioned techniques for thermal 
optimization with state-of-the-art techniques for finite-element based struc- 
tural optimization to produce a coupled thermal -mechanical optimization 
procedure. 

It is recognized that a complete and practical structural synthesis method 
should include design-oriented transient thermal analysis capability and the 
associated treatment of time-dependent constraints. Design variables for ther- 
mal control (insulation thickness, parameters of heat pipes and convective 
cooling systems) as well as parameters of thermal -stress relief devices 
(expansion joints) should be included in the procedure. Although consideration 
of such capability is discussed in the present paper, incorporation of these 
features into an optimization procedure is beyond the scope of the present work. 



In the present paper, a simplified thermal -mechanical optimization procedure is 
presented based on steady-state structural temperatures and design variables 
consisting of thicknesses of membrane elements and cross-sectional areas of bar 
elements. The procedure is based on a standard nonlinear mathematical program- 
ming formulation using the Sequence of Unconstrained Minimizations Techniques 
(SUMT) in which design requirements are incorporated with an exterior penalty 
function. Calculations of temperatures and stresses during resizing are per- 
formed by finite element analyses using one- and two-dimensional elements. 
Design requirements include limits on element stresses, temperatures and mini- 
mum gage. The procedure is demonstrated by the minimum-mass design of a wing 
box under combined mechanical loading and external heating. Designs are 
obtained with and without temperature constraints in order to indicate the 
effects of these requirements. 


SYMBOLS 


A 


mm 


(Bl 

tc x l, [C y l, 

E 


g 

H 


h 

h • 
mi 

[«3 

[K1 


k 

{L} 


area of a bar element 
bar minimum gage area 

row vector relating stress in a bar to displacement 

[Cl row vectors relating membrane stresses to displacements 
xy 

Young's modulus 

factors relating membrane thermal stresses to temperatures 
design constraint 

convective heat transfer coefficient 

membrane thickness 

membrane minimum gage thickness 

stiffness matrix 

conductivity matrix 

thermal conductivity 

structural load vector 
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nb 


nc 

nm 


P 

Q» q 
r 

{T> 

T 

{t} 

{U} 

{u} 

U, v, w 


V 


X-, Y-, Z- 
x, y 
a 


a 



T 


thermal load vector 

mass of structure 

forces per unit length (fig. 1) 

number of bar elements in finite element model 

number of design constraints 

number of membrane elements in finite element model 

pressure 

heat flux 

penalty function weighting factor 

vector of structural temperatures 

temperature of a bar or membrane element 

vector of design vari ables" 

displacement vector for structure 

displacement vector for bar or membrane element 

displacement components in X-, Y-, and Z-di recti ons, 

respectively 

x y x y xy J S Von Mises stress measure 
global coordinate directions (fig. 2) 

Cartesian coordinate directions in plane of membrane element 
coefficient of linear thermal expansion 
penalty function 
bar stress 

membrane stress components 
time 
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Subscripts: 

a allowable 

b bar 

1 index indicating discrete time or design variable (Appendix B) 

j index indicating constraint number 

m membrane 


DESIGN METHOD 
Overall Approach 

The objective of the design method is to determine a structural design 
having minimum mass and satisfying the requirements that the structural ele- 
ments have stress levels within specified 1 imits^ te mperati ires bel ow specified 
values, and structural sizes above minimum values. The approach is to formu- 
late and solve a nonlinear mathematical programming problem based on the 
Sequence of Unconstrained Minimizations Technique in which constraints are 
represented by an exterior penalty function (ref. 6). Calculations of stresses 
and temperatures are performed using finite element analyses. 

Structural Analysis 

Calculations of displacements and stresses are carried out using bar 
finite elements and triangular membrane elements having mid-side grid points 
(ref. 7). Pertinent equations used in the structural analysis and needed for 
gradient calculations are as follows: 

Equilibrium equation: 


Bar constitutive equation: 


[KJ {U> = {L} 


% = [B] {u) b ' E b a b T b 


( 1 ) 

( 2 ) 
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Membrane constitutive equation: 



The measure of stress in the membranes is the Von Mises stress V, given by 

v= [ 0 2 x + o 2 y -V y + 3aV* < 4 > 

Thermal Analysis 

Temperatures are calculated using the same finite element model used for 
the structural analysis. One-dimensional elements incorporate conduction along 
the length and also allow for convective heat transfer and constant heat loads 
applied to the element. The temperature is assumed to vary linearly along the 
element length. The two-dimensional elements incorporate conduction in the 
plane of the element and allow for heat loads normal to the plane of the element 
as well as normal to the sides of the element. The temperature is assumed to 
vary quadratically in the plane of the element and the temperature is uniform 
through the thickness. 

The equation for calculating structural temperatures is: 

[K] (T) = {L} (5) 

where [K] is the conductivity matrix for the structure, {T} is the vector 
of grid point temperatures, and {[} is the thermal force vector for the struc- 
ture including contributions from convective heating and prescribed heat flux. 
The conductivity matrices and thermal load vectors for membrane and bar elements 
are derived for a general set of thermal loads in Appendix A. 
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Optimization Method 


Problem definition and solution . - The optimization problem consists of 
determining a vector of design variables {t} which minimizes the mass of the 
structure m, subject to a set of constraints {g}. The problem is solved by 
the Sequence of Unconstrained Minimizations Technique in which constraints are 
incorporated by an exterior penalty function. Thus, the problem is formally 
posed as follows: 

nc 

Minimize $ s m + r ^ <g.> 2 (6) 

J v -1 

for a sequence of increasing values of the weighting factor r. The uncon- 
strained minimizations for each value of r are carried out using a first- 
order Quasi -Newton method denoted Davi don's Second Method (ref. 8). Similar 
to the method of Davidon-Fletcher-Powell (ref. 9), Davidon's Second Method 
requires only one evaluation of the function and gradient per iteration, and 
does not require a one-dimensional search to obtain the function minimum. The 
notation <g -> is defined as follows: 

J 



A 

(O 

C_j. 

V 

Jo if gj > o 

K 1f 9j <0 

(7) 

Constraints. - Each 

design variable (area of bar or thickness 

of membrane) 

is subjected to stress, 

temperature, and minimum gage constraints. 

These are 

summarized as follows: 




Bar stress: 

g = 1 

- ct/ct > 0 

3 

(8) 

Bar minimum gage: 

9 = 1 

- W A ^° 

(9) 
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I 


Bar temperature: 

g = 1 - T/T a > 0 

(10) 

Membrane stress: 

g = 1 - V/o & > 0 

(11) 

Membrane minimum gage: 

g = 1 - h. /h > 0 
3 min — 

(12) 

Membrane temperature: 

g = 1 - T/T a > 0 

(13) 


These functions are defined so that positive values correspond to satisfied 
constraints and negative values correspond to violated constraints. The total 
number of constraints is three times the number of structural elements. Expres- 
sions for components of the gradient of the constraints required during the 
optimization are derived in Appendix B. 

APPLICATION OF THE PROCEDURE TO MINI MUM MAS S DESIGN OF WING BOX 

Model Description 

The optimization procedure is demonstrated for a structure which requires 

finite element analysis but is simple enough to avoid a host of modeling 

details. The structure, loading, and design requirements have been selected 

to demonstrate salient features of the techniques and ideas developed in this 

paper. The overall structural configuration shown in figure 1 consists of 

upper and lower covers, 4 ribs and 3 spars. The loading includes ir, -plane 

forces per unit length N , N , and N , as well as pressure over the upper 

x y xy 

and lower surfaces. In addition, a uniform heat flux is applied over the 
lower surface of the structure. Boundary conditions include prescribed tem- 
peratures on edges 1 through 4 and prescribed zero displacements on edges 3 
and 4 (fig. 1). The thermal design requirements consist of an upper limit of 
1090 K (1500°F) on the temperature of all structural elements. The allowable 
stress for all structural elements is 0.875 GPa (127 ksi). Minimum gage 

A 

limitations are bar areas of 0.0064 cnr (0.01 in?) and membrane thicknesses of 
0.254 mm (0.01 in.). A summary of loads, boundary conditions, and material 
properties is given in table 1. 

The finite element model for the structure is shown in figure 2. The ribs 
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and spars are modeled by trusses and the covers are modeled by membranes. 
Modeling details are shown in figure 3. The model consists of 68 bars, 8 mem- 
branes and 30 grldpoints resulting in an optimization problem with 76 design 
variables. Coordinates of the grldpoints are listed in table 2 and the forces 
resolved Into grldpolnt forces are given in table 3. The boundary conditions 
expressed in terms of the gridpoint temperatures and displacements are given 
in table 4. 


Determination of Optimum Design 

Design with temperature constraints included . - The starting point for the 
optimization was a structure in which all structural elements were sized at 
minimum gage. This initial design had a mass of 39 kg (86 Ibm) and had a large 
number of constraint violations. The optimum design was obtained after 12 
unconstrained minimizations and had a mass of 486 kg (1070 Ibm). The cross- 
sectional areas of the bars and thicknesses of the membranes in the optimum 
design are listed in table 5. The thicknesses of the membranes in the upper 
and lower skins are shown in figure 4. The thickest element of tne skin 
(element 5) is in the lower surface adjacent to the built-in edge where the 
largest stresses occur, and this element is stress-critical. The next thickest 
element (element 7) is located at the interior of the lower surface where the 
highest temperatures occur, and this element is temperature-critical. A sunmary 
of the critical constraints is given in table 6(a). Of the 8 membrane elements 
used to represent the upper and lower skins, 3 are stress-critical, 3 are tem- 
perature-critical (all on the lower skin), and element 3 on the upper surface 
is sized at minimum gage. Element 2 has no active constraints. 

Design with temperature constraints omitted . - In order to assess the 
effect of including temperature constraints, the wing box calculations were 
repeated with temperature constraints (Knitted. The starting design was a 
minimum-gage structure. The final design was obtained after 11 unconstrained 
minimizations and had a mass of 280 kg (617 Ibm). This design is tabulated in 
table 7, and a list of the critical constraints is given in table 6(b). The 
distribution of skin thicknesses and identification of active constraints are 
shown in figure 5. The largest proportion of structural material is repre- 
sented by element 5 in the lower skin. The most notable difference in the skin 
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designs with and without temperature constraints is the size of element 7 which 
was temperati:re-critical when the upper limit of 1090 K (1500°F) was enforced. 
In the present design, the temperature of eltment 7 is 1133 K (1580°F) and the 
element is significantly thinner. Omittino the constraint permitted the higher 
temperatures and reduced the need for element 7 to conduct heat to cooler parts 
of the structure. The differences between the masses of the designs with and 
without temperature constraints is largely due to the di 'ference in mass of 
element 7. Membranes 6 and 8 which v,'ere also temperature-critical in the pre- 
vious design, had temperatures in the current design of 1109 K (1536°F) and 
1114 K (1545°F), respectively. 

CONSIDERATIONS FOR ADVANCED APPLICATIONS 

Advanced applications of the design procedure of this paper require 
development of the following features and capabilities: 

(1) design-oriented transient thermal analysis 

(2) treatment of time-dependent constraints 

(3) incorporation of effective temperature control devices 

(4) incorporation of thermal -stress relieving mechanisms. 

Each of these items is discussed in the succeeding paragraphs. 

Design-Oriented Transient Thermal Analysis 

A most pressing need before optimization procedures can be applied to 
practical aerospace structures under realistic load situations is the incor- 
poration of transient thermal analysis. While there are a variety of algorithms 
for tracing out time histories of temperatures (for example, ref. 10), such 
algorithms are typically iterative. Consequently, an iteration loop for ther- 
mal analysis must be embedded ‘n the optimization iteration loop leading to 
potentially long and expensive computer run<=. Design-oriented transient ther- 
mal analysis, perhaps based or Taylor series approximations in a manner similar 
to the analysis technique of references 1 and 11, has potential for efficient 
thermal optimization procedures. Additionally, the idea of updating tempera- 
tures only after selected design iterations (suggested in ref. 4) would 
reduce the total number of calculations. This technique should be especially 
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useful in the later stages of design where relatively small changes in the 
design variables occur and consequently, only small changes in the temperatures 
are expected. 


Time-Dependent Constraints 

Because temperatures and stresses can be time- dependent, behavior con- 
straints can be similarly time-dependent. Techniques are, therefore, needed 
for constraints such as 


g(x) - 1 - T(t)/T_ 


04) 


where t denotes time and T fl is the allowable temperature. One method of 
dealing with this type of constraint is to replace it by the integrated average 
s in reference 4. Thus, the actual constraint on time-dependent temperature 
is replaced by 



(15) 


where 



(16) 


and Tf is the length of time over which the temperature is being r >nitored. 
As pointed out in reference 4, one drawback to this approach is thrc it tends 
to smooth the penalty function and thereby causes temperatures to bo somewhat 
insensitive to changes in the design variables. An alternative approach is to 
satisfy the constraints at a finite number of discrete times (NT). Thus, the 
continuous time-varying constraint of equation (14) is replaced by the follow- 
ing constraints: 


g i (x i ) = 1 - T(x i )/T. 


(i = 1, 2, . . ., NT) 
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This approach has the disadvantage of Increasing the number of constraints that 
need to be satisfied. The advantage of this approach is that careful control 
of the choice of the discrete time values can be exercised. Thus, after a 
design based on an initial choice of discrete times is obtained, a check can 
be made to see if violations occur between the discrete time values. If there 
are violations, the design may be repeated with a refined choice of t. in 
the region of the violations. 

Incorporation of Temperature Control Devices 

Changing areas of bar elements and thicknesses of membrane elements is a 
relatively weak thermal control procedure. It is recognized that practical 
control of heated structures requires more effective the. mal control devices 
such as forced- convection active cooling systems (re f . 12). The number of 
possible design variables necessary to characterize tr»ese systems is potentially 
large and includes, for example, the sizes, spacing and shape of coolant pas- 
sages. The mass penalties include structural mass, coolant inventory mass as 
well as pumping penalties. Incorporation of this type of thermal control is 
desirable for effective thermal optimization, but appropriate characterizations 
of these systems are required before they can be incorporated into an overall 

optimization scheme for thermal /mechanical design. 

Thermal Stress Relief Mechanisms 

In reference 3, a method was developed for identifying situations where 
fully-stressed design procedures for structures under prescribed temperatures 
and mechanical loads would fail unless thermal stresses are reduced by means 
other than resizing. In that reference, a factor $ was computed as the 
largest fraction of the thermal stress that could be accommodated in each ele- 
ment. No attempt was made either to relate 3 to a physical thermal stress 
reduction mechanism or to assign a mass penalty to the value of $• 

Thermal stress reduction can often be achieved by allowing expansion at 
potentially highly thermally-stressed portions of structures. For example, 
expansion joints have been proposed in the design of thermally- loaded piping 
systems (ref. 13), and candidate thermal protection systems for space trans- 
portation-type vehicles are typically attached to the substructure by flexible 
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supports (ref. 14). Panels of tubular construction have been proposed for 
high- temperature applications because the curved shapes allow thermal expansion 
without significant thermal stresses (ref. 15). Relating 3 to an appropriate 
expansion device with an appropriate penalty could enhance the thermal /mechani- 
cal optimization procedure. 


CONCLUDING REMARKS 

This report describes a finite-element based methodology for minimum-mass 
design of structures subjected to combined thermal -mechanical loading and both 
strength and thermal requirements. To simplify the development, temperatures 
were assumed to be steady-state and controlled by structural sizes only. A 
mathematical programming method based on the Sequence of Unconstrained Minimi- 
z. ,A *ons Technique (SUMT) was used in which design requirements are represented 
by i exterior penalty function. Temperatures were treated as behavior vari- 
ables rather than fixed load-type quantities and temperatures were updated 
during resizing by finite element analyses. 

Design calculations were performed for a wing box with both mechanical 
loading and external heating and subject to design constraints on stress, 
minimum gage, and temperature. Optimum designs without temperature constraints 
were also obtained. Temperature constraints had a significant effect on both 
the distribution of structural material and the total mass in the final design. 

Some additional developments beyond the scope of the present work but 
needed for design of practical aerospace structures under realistic load situa- 
tions are identified. Among these developments are design-oriented transient 
thermal analysis capability, treatment of time-dependent constraints, and 
incorporation of effective temperature control devices and thermal -stress 
relieving mechanisms. 



APPENDIX A 


THERMAL ANALYSIS FINITE ELEMENTS 

The conductivity matrix and the thermal force vector for the bar and tri- 
angular finite elements are obtained by standard finite element methods (ref. 
16). It is necessary to start with a functional expressed in terms of tempera- 
tures which, when minimized according to the usual calculus of variations 
methods, yields the matrix equations governing the temperature distribution. 

The functional appropriate to three-dimensional heat conductions, including 
convection and surface heat-load effects, is 

u= ^l 7T ml 7T l dV+ /[ H (i- TT ») - qT ] dS <A1) 

where I k ] thermal conductivity matrix 

T temperature 

{ VT> temperature gradient vector = 

V volume of finite element 

H convective heat transfer coefficient 

T^ ambient temperature at surface of finite element 
S surface of finite element 

q heat flux over surface of element 

Bar Element 

The functional of equation (Al) for a one-dimensional element (bar) has 
the following form: 
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{ [i* (6)*" (?-"■) -*]*' 

V (y - V.i) ♦ V (r - Vi) 


- q-jATj - qgATg + H 1 


(A2) 


The terms in equation (A2) are defined with reference to figure 6(a). 

A cross-sectional area of bar 

P perimeter of bar cross-section 

To»* T w j, T^ ambient temperature adjacent to lateral surface of bar and end 
points 1 and 2, respectively 

Tj and Tg temperatures at bar end points 1 and 2, respectively 

H, H 1 , Hg convective heat transfer coefficient for medium adjacent to 

lateral surface of the bar and at end points 1 ar.d 2, 
respectively 

q» q-j* q 2 heat flux at lateral surface of bar and end points 1 and 2, 

respectively 


The temperature distribution along the bar is assumed to be linear, thus 



Substituting equation (A3) into 

(A2) gives 



m T iio b m - {T) T {L> b 

(A4) 

where . . 


m -R! 

(A5) 
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[K) b = 


M . hp& m . 

1 3 + H 1 A 


kA , HP*. 


l 6 


kA t HP*. 
l 6 


kA , HP*. . I* a 

T T" 2 


(A6) 


HP*T „0 

— + f + <h A + 


tL >b = 


HPU 


‘fn/t h 2 at 


°°2 


(A7) 


The matrix [K] b and the vector {L> b denote the conductivity matrix and 
the thermal force vector, respectively, for the bar element. 


Triangular Element 

The functional in equation (Al) for two-dimensional heat transfer over a 
triangular region has the following form: 
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where 



thermal conductivities with respect to coordinate 
directions in the plane of the element (fig. 6(b)) 

8T 

component of temperature gradient in x-direction 
ST 

— , component of temperature gradient in y-direction 
ay 


H, H.|, convective heat transfer coefficients corresponding to 

the media adjacent to the surface of the element and 
the three edges of the element, respectively 

Too* Tool* T^* To ,3 ambient tenperatures corresponding to the surface of 

the element and edges 1, 2, and 3, respectively 


*1* 9"] » ^2* ^3 


A 


heat flux normal to surface of the element and normal 
to edges 1, 2, and 3, respectively 

area of the triangular element 


S| » Sg, S3 edges of the triangular element 

Equation (A8) may be written more compactly as: 



(A9) 
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where 


1 

‘h 

0 

0 ’ 

[k] = 

0 

k x h 

V 


0 

k h 

kh 


m 

xy 

y 


(A10) 


q t HT 
n 00 

0 

0 (All) 

The temperature distribution is assumed to vary quadratically over the element 
as follows 

T = ag + a,x + a^y + a 3 xy + a^x 2 + a g y 2 (A12) 

The coefficients in (A12) may be expressed in terms of the temperatures and 
x, y coordinates at the six grid points of the element. This relationship is 
written in matrix form as 

(a) = [Q] {T> (A13) 

where 

{a} T = |a 0 a, a 2 a 3 a 4 a 5 J 
{T)T * [l T 2 T 3 T 4 T 5 T 6j 

I 1 x i x i*i x ? 

1 X 2 y 2 x 2 y 2 x| 

1 x 3 y 3 x 3 y 3 x| 

1 x 4 y* x*y 4 xl 

1 x 5 y s x 5 y 5 xf 

1 x 6 y 6 x 6 y 6 x| 
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Substitution of equation (A12) Into (A9) and use of equations (A13) through 
(A15) yields: 


U = i {T> T [K] m {T} - {T> T <L> m 


where 


[Kl = [Q ] J f tXl T [k] [X] dx dy 
si? j . 


[Q] 


)L} m = [Q l 1 / [X] T {s}dx dy 


'A 

+ h [Q] 


T + ds l + <<*2 + H 2 T {Y} T 

♦ ^3 + H 3 T »3>jf {Y}T ds 3^ 


ds. 


1X1 = 


x y 

1 0 

0 1 


xy 

y 

x 


x y 

2x 0 

0 2y 


(A17) 


+ {Y} T {Y} ds 1 + ^ J (Y> T {Y> ds 2 + ~~j (Y} T {Y> dSgj 

\ s 1 s 2 S 3 / 


[Q] 

(A18) 


(A19) 


(A20) 


{Y} T = [lxyxyx 2 y 2 J (A21 ) 

The Integration indicated by / is performed in closed form using tri- 
angular integration according to formulas given in reference 17. The line 
integrals indicated by /•/ , and J" are performed in closed form 


*1 *'s 2 


making use of the equations of the straight lines defining sides 1. 2, and 3. 
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The conductivity matrix and thermal force vector for the triangular ele 
ment are identified as iKl m and {L} m , respectively. Having derived the 
conductivity matrix and thermal force vector for the bars and membranes. It 
Is a simple matter to assemble the corresponding matrix and vector for the 
complete structure by the standard finite element method. The resulting 
matrix and vector are denoted C K 3 and {L} as used in equation (5) of the 
main text. 
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APPENDIX B 


GRADIENT OF PENALTY FUNCTION 


Expressions for the components of the gradient of the penalty function 
(eq. 6) are summarized In this appendix. The required gradient Is 


V <t> 


(■ 


V |m + r 


nc 


E 



(Bl) 


where m is the total mass of the structure and the gj are constraints 
defined in equations (8)-(13) of the main text. For convenience in notation* 
the vector of design variables {t> is defined as follows: 



{A} 

in} 


nb x 1 
nm x 1 


(B2) 


where {A} 
{h} 


is the vector of nb bar cross-sectional areas 
is the vector of nm membrane thicknesses. 


Thus, 


n 

m = m ({t}) = 2 Pi Jt i ti (B3) 

i=l 

where p. is the density of the i-th design variable and is the length 
of a bar element or planform area of a membrane element. Combining equations 
(Bl) and (B3) leads to the following expression for the i-th component of the 
gradient: 
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nc r jg i 

H. - Ol*1 ♦ 2r £ I Jtf | {< 9j >> 

1 J «1 *- 1 J 


(B4) 




The remaining task Is to obtain expressions for These expressions are 

summarized for each of the six types of constraints (eqs. (8)-(13)). 


Bar Stress Constraints 


Using equations (1), (2), and (8) 


M’M/ (?), ** 


(3 = 1,2 nb) (B5) 


where |ff"| Is the appropriate 6 x 1 sub-vector of 


|S|"«-'(|ft|-[K] <»>) 


associated with the two ends of the bar 


3T 

and ~ is the appropriate 2 x 1 sub-vector of 




associated with the temperatures at the two ends of the bar. 


(B6) 


(B7) 
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Bar Minimum-Gage Constraints 


Using equation (9): 

" a 

3t* " .2 °1.j-nb 

i t i 


(j = nb+1, nb+2, ...» 2nb) (B8) 


where 6^ Is the Kronecker delta. 


Bar Temperature Constraints 


Using equation (10): 


Hi. J.?ij-2nb 

3*1 - - T a at. 


(j = 2nb+l, 2nb+2 3nb) (B9) 


Membrane Stress Constraints 


Using equations (3), (4), and (11): 


= . _i!!.J-3nb (j = 3nb+l, 3nb+2, . 

9t- a 3t. 

1 a 1 

where, for pyery membrane element 


3V _ 
3t " 



+ 




3nb+nm) 


(BIO) 


(B11) 
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Membrane Temperature Constraints 

From equation (13): 


_ 1 *V3nb-2nm 

3t i ' " T i 3^ 

(j = 3nb+2nnH-l, 3nb+2nm+2, . . .» 3(nb+nm)) ( B14) 
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TABLE 1. - SUMMARY OF APPLIED LOADS, BOUNDARY 
CONDITIONS AND MATERIAL PROPERTIES FOR WING BOX 

(a) LOADS 


■ '■ ' 

UPPER 

SURFACE 

LOWER 

SURFACE 

N x 

N/m 

-38* 

*452 

Ibf/in 

-336 

-4000 





N y 

N/m 

-136 

147 

lbf/in 

-1200 

1300 





N Xy 

N/m 

60 

14 

lbf/in 

528 

128 





P 

Pa 

276 

6895 

psi 

.04 

1.0 





Q 

Watt/m 2 

0 

1.634 

BTU/in 2 sec 

0 

1 x 10- 6 


♦(minus sign compression) 
(b) BOUNDARY CONDITIONS 


Edge 1 

T = 987 K (1300°F) 

Edge 2 

T = 1061 K (1450°F) 

Edge 3 

T = 916 K (1270°F) , u = v = w = 0 

Edge 4 

T = 1047 K (1425°F) , u = v = w = 0 


(c) MATERIAL PROPERTIES 
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TABLE 2 


GRID POINT COORDINATES OF WING BOX 



1 















TABLE 3 


APPLIED MECHANICAL FORCES AT GRID POINTS 


OF WING BOX FINITE ELEMENT MODEL 


Grid 

Px 

p y 

P 

Z 

Point 

N 

lbf 

N 

lbf 

N 

lbf 

1 

- 11387 

- 2560 

39856 

8960 

- 49 

-11 

2 

-125262 

- 28160 

248686 

64000 

- 71 

-16 

3 

- 51244 

- 11520 

102487 

23040 

- 22 

- 5 

4 

-244830 

- 55040 

- 69499 

-15624 

592 

133 

5 

- 30372 

- 6828 

-308404 

-69332 

1779 

400 

6 

229635 

51624 

- 84694 

-19040 

1188 

267 

7 

79712 

17920 

-125262 

-28160 

- 71 

-16 

8 

0 

0 

0 

0 

-142 

-32 

9 

- 79712 

- 17920 

125262 

28160 

- 71 

-16 

10 

-948948 

-213332 

30372 

6828 

1779 

400 

11 

0 

0 

0 

0 

3559 


12 

948948 

213332 

- 30372 

- 6828 

1779 

400 

13 

39856 

8960 

- 62631 

-14080 

- 71 

-16 

14 

0 

0 

0 

0 

-142 

-32 

15 

- 39856 

- 8960 

62631 

14080 

- 71 

-16 

16 

-474483 

-106668 

15177 

3412 

1779 

400 

17 

0 

0 

0 

0 

3559 

800 

18 

474483 

106668 

- 15177 

- 3412 

1779 

400 

19 

79712 

17920 

-125262 

-28160 

- 71 

-16 

20 

0 

0 

0 

0 

-142 

-32 

21 

- 79712 

- 17920 

125262 

28160 

- 71 

-16 

22 

-948948 

-213332 

30372 

6828 

1779 

400 

23 

0 

0 

0 

0 

3559 

800 

24 

948948 

213332 

- 30372 

- 6828 

1779 

400 


29 















TABLE 4. - MECHANICAL AND THERMAL BOUNDARY 



CONDITIONS FOR WING BOX MODEL 


3i 


i] 










TABLE 5. - FINAL DESIGN OF WING BOX INCLUDING 
STRENGTH AND TEMPERATURE CONSTRAINTS 
(a) BAR AREAS 


Area ] 

Bar 

Are< 

l 

Bar 

Area 

cm 2 

in 2 

cm 2 

in 2 

cm 2 

in 2 

.6232 

.0966 

24 

.1639 

.0254 

47 

.1135 

.0176 

.1284 

.0199 

25 

.6200 

.0961 

48 

.1439 

.0223 

.5871 

.0910 

26 

.0684 

.0106 

49 

.0645 

.0100 

.1548 

.0240 

27 

.3684 

.0571 

50 

.0877 

.0136 

.1381 

.0214 

28 

1.7123 

.2654 

51 

.4994 

.0774 

.0916 

.0142 

29 

.3568 

.0553 

52 

.0710 

.0110 

.1574 

.0244 

30 

1.7445 

.2704 

53 

.4471 

.0693 

.0690 

.0107 

31 

2.2090 

.3424 

54 

.2239 

.0347 

2.2529 

.3492 

32 

.1194 

.0185 

55 

.0652 

.0101 

.3619 

.0561 

33 

.1303 

.0202 

56 

.1084 

.0168 

.1374 

.0213 

34 

.0806 

.0125 

57 

.8426 

.1306 

2.6697 

.4138 

35 

.5174 

.0802 

58 

.0942 

.0146 

.0794 

.0123 

36 

.3948 

.0612 

59 

.0665 

.0103 

.0742 

.0115 

37 

.1832 

.0284 

60 

.0742 

.0115 

.0652 

.0101 

38 

.1529 

.0237 

61 

.0877 

.0136 

.0755 

.0117 

39 

.1013 

.0157 

62 

1.8413 

.2854 

.1161 

.0180 

40 

.0826 

.0128 

63 

.1032 

.0160 

.0774 

.0120 

41 

.0652 

.0101 

64 

.1174 

.0182 

.8600 

.1333 

42 

.1123 

.0174 

65 

.4555 

.0706 

.2748 

.0426 

43 

.0852 

.0132 

66 

.3490 

.0541 

.1974 

.0306 

44 

.1768 

.0274 

67 

2.4032 

.3725 

.0819 

.0127 

45 

.0774 

.0120 

68 

2.6858 

.4163 

.4568 



.0708 

46 

.3703 

.0574 





(b MEMBRANE THICKNESSES 


Membrane 


Thickness 



.0101 
.0269 I .0106 


Membrane 



Thickness 

cm 

in 

1.3861 

.2769 

.7488 

.0935 

.5457 

.1090 

.2948 

.0368 


Final Mass - 485 kg (1070 lbm) 
































TABLE 6. - STATUS OF CONSTRAINTS IN WING BOX DESIGN 


(a) WING BOX WITH TEMPERATURE CONSTRAINTS 

Stress C 

Critical 

Temperature Critical 

Bars 

Membranes 

Bars 

Membranes 

19 

1, 4, 5 

19, 20, 23, 
24, 27, 

33, 34, 35, 
42, 43, 

49, 50, 51 

6, 7, 8 


(b) WING BOX WITHOUT TEMPERATURE CONSTRAINTS 

Stress 

Critical 

Temperature Critical 

Bars 

Membranes 

Bars 

Membranes 

19, 27, 63, 

2, 5, 6, 8 

N.A. 

N.A. 

64 


























TABLE 7. - FINAL DESIGN OF MING BOX 
WITHOUT TEMPERATURE CONSTRAINTS 

(a) BAR AREAS 



.0645 

.1742 

.0645 

.0774 

.1252 

.2361 

.0677 

.0748 

.2574 

.0910 

2.883 

3.974 

.2858 

.3884 

.0671 


.0120 

.0194 

.0366 

.0105 

.0116 

.0399 

.0141 

.4472 

.6165 

.0443 

.0602 

.0104 



Area 


Area 

Bar 

cm 2 

in 2 

Bar 

cm 2 

in 2 

24 

.2935 

.0455 

47 

.0729 

.0113 

25 

.1981 

.0307 

48 

.0658 

.0102 

26 

.2013 

.0312 

49 

.0735 

.0114 

27 

1.548 

.2398 

50 

.0652 

.0101 

28 

.0729 

.0113 

51 

.0664 

.0103 

29 

.0910 

.0141 

52 

.0645 

.0100 

30 

.3032 

.0470 

53 

.3439 

.0533 

31 

.0652 

.0101 

54 

.0677 

.0105 

32 

.0748 

.0116 

55 

.2310 

.0358 

33 

.0761 

.0118 

56 

.4942 

.0766 

34 

.2052 

.0318 

57 

.0845 

.0131 

35 

.2806 

.0435 

58 

.0684 

.0106 

36 

.2052 

.0318 

59 

.1013 

.0157 

37 

.0645 

.0100 

60 

.0658 

.0102 

38 

.2858 

.0443 

61 

.0710 

.0110 

39 

.2632 

.0408 

62 

.1955 

.0303 

40 

.0852 

.0132 

63 

.1864 

.0289 

41 

.0671 

.0104 

64 

.2426 

.0376 

42 

.0710 

.0110 

65 

.4968 

.0770 

43 

.0910 

.0141 

66 

.0658 

.0102 

44 

.1974 

.0306 

67 

.1413 

.0219 

45 

46 

.0664 

.0645 

.0103 

.0100 

68 

.1600 

.0248 


b) MEMBRANE THICKNESSES 


Membrane 


Thickness 


Membrane 


Thickness 



290 

395 

387 


.0331 


Final Mass = 280 kg (617 lbm) 

































Loads on Upper Surface 


P 



Loads on Lower Surface 

Figure 1.- Wing box model showing applied loads and edges where 
boundary conditions are specified (see table 1). 
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Spar 1 Spar 2 Spar 3 


Figure 2 Finite element mode! of wing box (for clarity not all 
elements are shown). 



















(a) Upper skin 

S indicates stress-critical 
T Indicates temperature-critical 
M indicates minimum gage 


in. cm 


0.6 

0 



Figure 4.- Distribution of membrane thickness in final design 
of wing box with temperature constraints included. 



in. cm 



S indicates stress-critical 
M indicates minimum gage 


in. cm 



<b) Lower skin 
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Figure 5.- Distribution of membrane thickness in final design 
of wing box without temperature constraints. 







H 2 ^2 

End 2 


(a) Bar Element 


3 



Side 1 

(b) Triangular Element 
Figure 6.- Thermal finite elements. 
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